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Abstract 



The expected root-mean-square value of a matrix element A a p in a 
classically chaotic system, where A is a smooth, ^-independent function 
of the coordinates and momenta, and a and (3 label different energy 
eigenstates, has been evaluated in the literature in two different ways: 
by treating the energy eigenfunctions as gaussian random variables and 
averaging |^4 Q /3| 2 over them; and by relating |^4 a ,a| 2 to the classical time- 
correlation function of A. We show that these two methods give the 
same answer only if Berry's formula for the spatial correlations in the 
energy eigenfunctions (which is based on a microcanonical density in 
phase space) is modified at large separations in a manner which we 
previously proposed. 
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Hamiltonian systems which are classically chaotic have quantum energy eigenvalues, 
eigenfunctions, and transition matrix elements which can be profitably analyzed statisti- 
cally M$M- Our focus in this paper will be on matrix elements (in the energy-eigenstate 
basis) of operators whose Weyl symbols are smooth, ^-independent functions of the classical 
coordinates and momenta. Two different methods have been proposed in the literature for 
calculating the root-mean-square statistical average of these matrix elements in the limit of 
small h. One method is to compute this average by treating the energy eigenfunctions as 
gaussian random variables; the other relates the average to the operator's classical power 
spectrum. Our goal is to see whether or not these two methods give the same result, a ques- 
tion which was first raised by Austin and Wilkinson ||. We find that the methods do agree, 
but only if our recently proposed modification of Berry's formula @] for the spatial correla- 
tions in energy eigenfunctions of chaotic systems is invoked when the spatial separation is 
large compared to any relevant classical distance scales in the problem ||. 

We begin by reviewing the power-spectrum method, essentially following the original 
arguments of Feingold and Peres ||; more rigorous treatments leading to the same result 
have been given by Wilkinson [0] and Prosen [Q . To simplify the discussion, we will consider 
hermitian operators which are functions of only the coordinates q (and not the momenta 
p). Given a suitable operator A of this type, we begin by defining 



dte- t /2T 'e iuJt (a\A t A\a) , (1) 

where A t = e ' lHt / n j\ e - tHt / n [ s the relevant operator at time t in the Heisenberg-picture, |a) 
is an energy eigenstate with energy E a , uj is a parameter, and r c is a time cutoff which may 
be needed for convergence of the integral. 

We now evaluate F in two different ways. First, we use Shnirelman's theorem |9|-|1"3 



which says that, in the limit of small h, the expectation value of an operator O in an energy 
eigenstate is equal to its classical, microcanonical average at the corresponding energy, 

(a|0|a) = J dn Ea O w (p, q) , (2) 

where 0\y{p, q) is the Weyl symbol of the operator O, and dfiE denotes the Liouville measure 
on the surface in phase space with energy E, 

1 d^p d^Q 

d » E = W)Wy 5{E ~ HwMy (3) 

Here / is the number of degrees of freedom, Hw(p, q) is the Weyl symbol of the hamiltonian 
operator H, and p(E) is the semiclassical density of states, 

m = /f^^-^p^))- ( 4 ) 

Note that d\iE is a purely classical object; the factors of h cancel between eqs. @ and (f|). 
Also, Shnirelman's theorem is proved for principal symbols instead of Weyl symbols, but 
there is no difference in the h — > limit which the theorem also requires; see |14| for a 
thorough discussion. 
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We now apply eq. (0) in eq. (|l|), making the approximation (valid in the h — > limit) 
that 

J dfi E (A t A) w = J dp E A w (q t )A w (q) , (5) 

where A w (q) is the Weyl symbol of the operator A (which, by assumption, depends only on 
q and not p), and q t is the classical coordinate at time t, assuming an initial point (p, q) 
on the surface with energy E in phase space. We therefore obtain 

F = [ +0 ° dt e"' 2 / 2 ^ / dfi Ea A w (qt) A w (q) . (6) 

J —oo J 

We now evaluate F in a different manner: we insert a complete set of energy eigenstates 
to get 

/+oo 9j , 

dte' t/2r 'e lult (a\A t \(3)((3\A\a) 

= E [ +0 ° dte- t2 / 2r *e l ^- E P +h ^ h A a pAp a 

= 27rnJ28 n/Tc (E a - Ep + fko)\A aP \ 2 , (7) 

where A a p = (a\A\f3), and 5 e (E) denotes a smeared delta function with a width of e. We now 
assume that each eigenstate has a random character, so that, with E a and Ep each varied 
over a small range, there is a smooth distribution of values for |y4 a/ 3| 2 . Let this distribution 
be characterized by an expected value which we will call (I^Iq,^ | 2 ) . If we also take the width 
H/r c of the smeared delta functions to be somewhat larger than the mean level spacing, 
equal to l/p(E) in the limit of small h, we can replace the sum over (3 in eq. (ffl) by an 
integral over p(Ep)dEp. Thus we have 

("OO 

F = 2nh / dE p p{Ep)8 1/p {E a - E p + hu)(\A a/3 \ 2 ) 

J 

= 27ihp{E a + huj)(\A a p\ 2 ) . (8) 

Equating the right-hands sides of eqs. ([[[) and (§) gives us the desired formula for (|A Q( g| 2 ); 
however, its accuracy to subleading order in h can be improved by symmetrizing on E a and 
E § to get 

(\A a(3 \ 2 ) = - / dte-^rte™ / dp E A w (q t )Aw(<l) , (9) 

Th J-oo J 

where E = \{E a + Ep) is the mean energy, %uj = Ep — E a is the energy difference, and 
th = 2-n;Tip(E) is the Heisenberg time. If we hold E and uj fixed in the limit of small H, the 
right-hand side of eq. (§) is simply 1/th times the classical power spectrum of the observable 
A at energy E, with any structure on frequency scales less than 27r/r H smeared out by the 
time cutoff. Eq. (|9|) is the first of two formulas for (IAq^I 2 ) which can be found in the 
literature. 
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We get the second formula [15]-|l7| by first writing the squared matrix element in terms 
of the eigenfunctions as 

\A a p\ 2 = J dfq'rM)A w (q!)M4) J dfqP p {q)A w (q)1> a (q) . (10) 

In a chaotic system, the individual eigenfunctions can be treated as independent random 
variables with a gaussian probability distribution PJfS|-2I]. Because it is gaussian, this 
distribution is completely specified by the two-point correlation function 

C(q',q|E) = (^(q')^(q)), (11) 

where the angle brackets denote averaging over the probability distribution for ip(q) given 
the energy E. Averaging eq. (|10"D over this probability distribution yields 



(I^VD = / d f (/d f qC{<i,q[\E a )A w (q[)C{q[,q\Ep)A w (q) . (12) 

This is the second formula for (|/l a/ 3| 2 ) which can be found in the literature. The question 
is whether or not it is the same as the first formula, eq. ([|). 

Of course, in order to answer this question we need an explicit expression for C(q, q.'\E). 
Berry [|J conjectured that, in the small-fr limit, 

C(cj,q\E) = J ^ e ^'-^(Q-I(q' + q )) , (13) 

where the Liouville integral is over (P, Q). However ||, this formula for C(q', q\E) appears 
to be too simple to be able reproduce the classical power spectrum of A which appears in 
eq. (§• 

In a separate paper ||, we have argued that, in fact, eq. (p~3|) must be modified whenever 
the separation |q' — q| is large, in the sense that the shortest classical path with energy E 
which connects q to q' is not well approximated by a linear function of time. This will be 
generically true in eq. (|9]), since both q and q' are integrated, and since the factors of A^y(q) 
and Ay/{c{) do not force q and q' to be close together. When q and q' are far apart, eq. ([13]) 
should be replaced with 

C ttAE) = -, Wo 2 fc w f+1 w 2 E \ D P \ 1/2 ™s[S p /h - (2u p + / - 1)tt/4] , (14) 



p(E)(2nh)U+ 1 )/ 2 



paths 



where the sum is over all classical paths connecting q to q' with energy E; each path has 
action S p = p-<iq, focal point number u p , and fluctuation determinant 



D p = det 



( d 2 s p d 2 s p 

9q9q' dEdq! 

d 2 S p d 2 S p 

\ dqdE dE 2 



(15) 



Eq. (0) actually holds only if the system is invariant under time reversal; otherwise a more 
cumbersome formula is needed ||. The final formula for (|y4 Q/ g| 2 ) turns out to be the same in 
either case, and so to simplify the notation we will use eq. (|i~4"D. Eq. (|i~4"|), or its replacement 
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for a system which is not time- reversal invariant, is valid as long as the contributing path 
of least action has S p /% ^> 1. This is of course true generically in the limit of small h. 

We now show that if we use eq. (fH| ) for C(q / ,q|£'), eq. ( p~2"D gives the same result for 
(\A af3 \ 2 ) as eq. ©. 

We begin by substituting eq. (|T4| ) into eq. (|12|). Since we are interested in the limit of 
small h with E and u held fixed, we can usually replace E a and Ep with E. We then make 
use of the "diagonal approximation" in which the double sum over paths is collapsed to 



a single sum. In related calculations |2^JT^], this can be justified by the rapidly oscillating 
phases of the off-diagonal terms as long as the single sum includes only those paths whose 
elapsed times are less than the Heisenberg time. We assume the same condition holds here. 
The product of cosines in each diagonal term then yields a single cosine which is slowly 
oscillating, and we get 



(|A^| 2 > = KE y {27ch y + i J dt</dfqA w (4)A w (q) £ \D p \ cos(ur p ) . (16) 
The sum is over paths from q to q' with energy E, and elapsed time 

(17) 



dS p 



p ~ dE 



E=E 



less than the Heisenberg time Th- We have implicitly assumed that the topological quantity 
u p does not change as the energy of the path is varied from E — \hio to E + ^Tiu. 
To make further progress we need to rewrite the fluctuation determinant as 



dp 


dr 


<9q' 


<9q' 


dp 


dr 


dE 


dE 



Z) p = det - 2 . (18) 

V dE dE I 

Here p = —dS v /d<\ is the momentum at the beginning of the path, and r = r p is the elapsed 
time along the path, given by eq. (0). With these definitions, eq. (|18D follows immediately 



from eq. (TT5|). Eq. (|Tj|) shows us that \D p \ can be thought of as a jacobian for a change 
of variables from the final position q' and total energy E to the initial momentum p and 
elapsed time r [pq| . To make use of this, we insert 1 = / dE 5{E — Hyy(p,q)) into eq. flTBD 



and change variables. We now have 

rTH r d-fpdfq 



Aa ^ 2) = ^MeY r dr I j£w ? 6{E ~ Hwip ' q)) ^(^^(^^(q) • (19) 



paths 



The sum is over all paths which begin at (p, q) and have elapsed time r. However, there is 
only one such path, and so the sum over paths may be dropped. Also, q' is the position at 
time r, and it is now more properly denoted q T . Using eq. (||) and r H = 2nTip(E), we see 
that eq. (|l9| ) can be rewritten as 

(lA^I 2 ) = — / drcos(^r) / dfJ, E A w (q T )A w (q) . (20) 
th jo j 

Using the fact that time-translation invariance implies that / efytg A\y{(\ t )A\y{c\) is an even 
function of r (even if the system is not time-reversal invariant), we see immediately that 
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eq. (pCf) is equivalent to eq. (§), up to the issue of the detailed treatment of the large-time 
cutoff. This is our main result. 

Another quantity of interest is the size of the fluctuations in the diagonal matrix elements 
A aa . If we first shift A (if necessary) so that (A aa ) = J dfiE a A w (q) = 0, then the object 
we wish to evaluate is (|A aQ | 2 ). This has been done previously [F2^],^,|nj by making use of 
the trace formula [^6j^9| , ^0| and properties of periodic orbits. Here we will compute (lAo^l 2 ) 
by averaging over the gaussian probability distribution for energy eigenfunctions [|H)|-|i"T|] . In 
the case of a system which is not invariant under time reversal, the energy eigenfunctions 
are generically complex, and the relevant formula is []23| 

WxWM = (»2><^4> + (WaWM , (2i) 

where ip^ = ip a ((ii). If the system is invariant under time reversal, the energy eigenfunctions 
are real, and we have instead [2H| 



Combined with the previous results for (|74 a/ g| 2 ), we find that 



9_ 



m 
-th 



dr J dfJ, Ea A w (q t )A w (q L ) . 



(22) 



(23) 



Here g = 2 for a system which is invariant under time reversal, and g = 1 for a system which 
is not. Eq. (p3|) is in agreement with the results of p7|j28|JT6[ . 
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